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ABSTRACT 

We present Monte Carlo simulations of the Diffuse Ha Galactic Background. Our 
models comprise direct and multiply scattered Ho; radiation from the kpc scaleheight 
Warm Ionized Medium and midplane H II regions. The scattering is off dust that is 
assumed to be well mixed with the gas, with an axisymmetric density distribution taken 
from the literature. The results of our simulations are all-sky Ha images that enable us 
to separate out the contributions of direct and scattered radiation. We also determine 
how far the model Ha photons have traveled, i.e., how far we see into the Galaxy at 
Ha. Our models reproduce the overall characteristics of the observed Ha background 
and predict the scattered Ha intensity at high latitudes is in the range 5% to 20% of the 
total intensity, in agreement with estimations based on [S II] /Ha and [O III] /Ha line 
ratio measurements. The polarization arising from dust scattering of Ha from midplane 
H II regions is predicted to be less than 1% at high latitude. 

Subject headings: diffuse radiation — H II regions — ISM — radiation transfer — dust 
scattering — polarization 
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1. Introduction 

Ubiquitous but faint Ha proves the presence of the 
Warm Ionized Medium (WIM), a major but poorly 
understood component of the ISM containing the ma- 
jority of the ionized gas in the Galaxy (see a review 
by Reynolds et al. 1995). The ionized gas is seen at 
all galactic latitudes and has a scale height of about 
900pc (judged from the dispersion measure of the ra- 
dio signals from pulsars in globular clusters), a tem- 
perature of around 10000K, and volume filling factor 
in the range 0.1 to 0.4. In the solar neighborhood, 
the total power represented by the ionizations pro- 
ducing the Ha is about 15% of the total ionizing rate 
from OB stars, or about 100% of the power ejected by 
supernovae. The source of ionization is not well un- 
derstood and has been narrowed down to either (a) 
the leakage of the ionizing photons from OB stars to 
the diffuse ionized gas, which actually provides the 
correct spectrum of the O, N, and S ions (Domgorgen 
and Mathis 1994; see also Dove & Shull 1994; Dove, 
Shull, & Ferrara 1999; Miller & Cox 1993; Elmegrecn 
1998); or (b) some unknown, more distributed source 
of ionization (Sciama 1995). 

Low spatial resolution mapping of the Ha emission 
(Reynolds 1990) has revealed the following character- 
istics. 

• At low latitudes (\b\ < 10°) the Galactic mid- 
plane is dominated by Ha emission from bright 
H II regions. 

• At intermediate Galactic latitudes (10° < \b\ < 
50°) the Ha intensity follows a esc \b\ law with 
the average intensity being < 7 a sin|6| >« 
1R, where R = 10 6 /47r Ha photons s _1 cm~ 2 
steradian -1 . 

• Above |6| = 50° the emissivity appears to be 
about a factor of two below the esc |6| law for 
lower latitudes. 

• Localized departures from the esc |6| law are ev- 
ident at all latitudes, indicating the clumpy and 
filamentary nature of the WIM. These features 
are becoming more evident through early re- 
sults from the Wisconsin Ha Mapper (WHAM) 
all sky survey (Reynolds et al. 1998; Haffncr, 
Reynolds, & Tufte 1998). 

The Galactic contributions to the Ha background 
are direct and scattered Ha from the WIM, direct Ha 



from midplanc H II regions, and Ha from H II regions 
that is scattered by dust at high latitudes. Observa- 
tions of [S II] /Ha line ratios along sightlines at low 
latitude (towards midplane H II regions) and at high 
latitudes in the diffuse ISM (Reynolds 1988), show 
that the line ratios in the diffuse ISM are larger by 
up to a factor of ten than those in traditional H II re- 
gions. This indicates that the ionization state of the 
WIM is very different from that of H II regions, and 
based on these observations, Reynolds (1990) con- 
cluded that the scattered component from midplanc 
H II regions was at most 15% of the total Ha inten- 
sity. 

Ferrara et al. (1996) presented a radiation transfer 
model for Ha emission and scattering in the edge-on 
galaxy NGC 891. They found that at 600pc from 
the midplanc, the scattered light from H II regions 
is around 10% of the total Ha intensity. Their mod- 
els also predict a dust scattered polarization of less 
than 1%. The analysis that we present in this paper 
is very similar to that of Ferrara et al. (1996) except 
we are now viewing the radiation transfer simulation 
from inside the Galaxy. Our goal is to determine, 
from a theoretical/modeling perspective, the relative 
contributions of direct and scattered Ha at high lat- 
itudes, how large a region of the ISM contributes to 
the observed intensity, and whether our models can 
reproduce the general characteristics of the Ha back- 
ground. In addition, we present a prediction of the 
polarization arising from dust scattering of the Ha 
emission. The ingredients of our models are presented 
in §2, and the scattered light models and polarization 
predictions are presented in §3 and §4. 

2. Models 

We perform the radiation transfer with a 3D Monte 
Carlo scattering code. We have previously modeled 
the polarization patterns and diffuse scattered light in 
external galaxies (Wood 1997, Wood & Jones 1997). 
These investigations adopted smooth axisymmetric 
models for the dust and illuminating starlight. In 
order to determine the propagation of starlight in 
clumpy media and simulate all-sky maps viewed from 
within a model galaxy, we have modified our codes to 
include the effects of scattering, clumping, dust scale- 
heights, and source distributions in a three dimen- 
sional galactic model. The ingredients of our mod- 
els are the Ha emissivity, dust density distribution 
and scattering parameters, and radiation transfer al- 
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gorithm, described below. 

2.1. Density Grid 

Our model galaxy density is constructed on a 3D 
linear Cartesian grid comprising 201 x 201 x 201 grid 
cells. The physical extent of our model galaxy is 
±15kpc, so each grid cell is a cube of w 150pc on 
a side. 

The number densities (cm -3 ) for the five phases 
of the ISM (molecular, cold neutral, warm neutral, 
warm ionized, and hot ionized) are determined from 
analytic expressions presented by Ferriere (1998, Eqs. 1 
through 9). This is a smooth axisymmetric model of 
the ISM. In the model at the solar radius the WIM has 
a scaleheight of lkpc and midplane number density 
of 0.0237cm -3 . The neutral hydrogen is composed 
of three components: Gaussian scaleheights of 127pc 
and 318pc and an exponential scaleheight of 403pc 
with the total midplane number density (cold neutral 
plus warm neutral) being 0.566cm -3 . We have also 
included the spiral arm density determined by Tay- 
lor & Cordes (1993). The spiral arms were inserted 
into our density grid using the FORTRAN subrou- 
tines that are described in the Appendix of Taylor & 
Cordes (1993). To complete our model ISM density 
we have inserted an evacuated region of 200pc radius 
centered on the Sun's location. This approximately 
represents the low density "local bubble" (e.g., Cox 
& Reynolds 1987), within the constraints of the grid 
size. 

We assume that the dust and gas are well mixed 
in all phases of the ISM and that the dust-to-gas ra- 
tio is constant throughout the Galaxy. We have not 
attempted to model any vertical or horizontal depen- 
dence of the dust to gas ratio within the Galaxy as 
was speculated by Ferrara et al. (1996). This would 
involve a detailed modeling of many datasets includ- 
ing IRAS, Ha, and H I and is beyond the scope of 
our current investigation to estimate scattered light 
levels. The dust scattering parameters are taken to 
be those of a standard ISM mixture (see §2.3). 

2.2. Ha Emissivity 

The Ha emissivity comprises emission from the 
WIM and that from H II regions. The WIM Ha emis- 
sivity is proportional to the square of the WIM den- 
sity (Eq. 7 in Ferriere 1998). In the neighborhood of 
the Sun and near the Galactic midplane, the mean 
density in the WIM is 0.025 cm -3 and the emissiv- 



ity is 6.7 x 10 -16 Ha photons/cm 3 /s (as derived in 
Reynolds 1990). Extrapolating this emissivity using 
the Ferriere density results in a total WIM emissivity 
from our model Galaxy of 10 52 Ha photons/s. 

Bright H II regions are the dominant source of Ha 
near the Galactic equator, with angular sizes ranging 
from arcminutes to tens of degrees. In order to ap- 
proximate the amount and distribution of Ha from 
local H II regions, we have utilized the Garmany et 
al. (1982) catalog of O stars within 2.5kpc of the Sun. 
This catalog (obtained electronically from the NASA 
Astrophysical Data Center) provides the spectral type 
of the sources and we convert this to a number of ion- 
izing photons from each source using Table 1 of Miller 
& Cox (1993). We assume that 6% of the ionizing 
photons from each source escape the Galaxy (Bland- 
Hawthorn & Maloney 1999) and 15% go into produc- 
ing the WIM (Reynolds 1990). Our escape fraction 
is consistent with the values of 6% and 3% (depend- 
ing on star formation history) estimated from models 
of Dove et al. 1999. The remaining 79% from each 
source are assumed to produce local "point source" 
HII regions with the Ha photon luminosity equal to 
0.46 x 0.79 Arizing = 0.36A ion i zing , assuming Case B 
recombination (Martin 1988). We have assumed that 
Nho — 0.36A ion i zin g for each source irrespective of 
location. This is clearly a simplification, since the 
Ha luminosity will depend on the local density which 
in turn determines what fraction of ionizing photons 
escape the Galaxy, what fraction contributes to pro- 
ducing the WIM, and what fraction is converted to 
Ha within the local H II region. This calculation is 
beyond the scope of our current investigation and we 
therefore globally adhere to the relative fractions we 
have quoted. The model of Miller & Cox (1993) at- 
tributed the WIM Ha emission to extended "Strom- 
gren volumes" created by OB associations and such a 
scenario would also include Ha emission from super- 
shells that may have formed around OB associations. 
Note that in our model we have separated the Ha into 
two components: a smooth component representing 
the extended emission of Miller & Cox and the point 
source H II regions with Ha luminosities as described 
above. 

For regions beyond 2.5kpc (where the Garmany et 
al. catalog is incomplete) we extend the distribution 
of point sources by placing sources randomly in the 
molecular ring (from Ferriere 1998) and in the spi- 
ral arms of Taylor & Cordes (1993). The (x,y) lo- 
cations of all sources in our simulation arc shown in 



3 



Fig. 1. We set the Ha luminosity of each of the addi- 
tional randomly placed point sources to be 10 49 Ha 
photons/s, (comparable to that of the Orion Nebula) 
and add enough sources such that the total luminos- 
ity of the H II regions in our simulation is equal to 
that of the WIM, L H n = L wlM = 10 52 Ha pho- 
tons/s. Several authors have estimated the relative 
contributions of Ha from the WIM and H II regions 
in nearby spiral galaxies, with the ratio found to be 
in the range 0.47 < Ha(H II)/Ha(TOTAL) < 0.7 for 
each galaxy as a whole (Ferguson, Wyse, & Gallagher 
1996; Walterbos & Braun 1994; Veilleux, Cccile, & 
Bland-Hawthorn 1995). The relative Ha luminosities 
for the H II regions and the WIM in our simulation 
have been chosen to match this observational con- 
straint. 



tion, HG(6), with g = 0.44 

HG{6) = i-(l - g 2 )/(l + g 2 - 2gcos9)^ 2 



(1) 




Fig. 1. — The x and y locations of the "point source 
H II regions." Locations of O stars from the Gar- 
many et al (1982) catalog and random point sources 
distributed in the molecular ring (Ferriere 1998) and 
in the spiral arms of Taylor & Cordes (1993). The 
Sun is located at (0, 8.5, 0). 



2.3. Opacity and Scattering Parameters 

We assume that the dust plus gas is represented 
by a Mathis, Rumpl, & Nordsieck (1977) mixture. 
The adopted parameters at Ha are a total opacity of 
220cm 2 g _1 , albedo a = 0.5, and a scattering phase 
function given by the Heyney-Greenstein phase func- 



The peak linear polarization, pi = 0.5, is taken from 
White (1979). For the adopted densities and opaci- 
ties, the total optical depth diametrically across the 
galactic equator is r cq « 50, while at the Sun's radius, 
the total optical depth perpendicular to the Galactic 
disk is 0.4. 

There has been considerable interest in determin- 
ing dust parameters in the diffuse ISM (Murthy & 
Henry 1995; Hurwitz 1994; Gordon ct al. 1994), 
reflection nebulae (Calzetti, et al. 1995), and pro- 
tostellar clouds (Whitney, Kenyon, & Gomez 1998). 
The effects of changing the albedo and phase func- 
tion asymmetry on scattered light patterns have been 
presented by other groups. Increasing or decreasing 
the albedo obviously increases or decreases the scat- 
tered light component. In their model of the UV dif- 
fuse Galactic light, Murthy & Henry (1995) showed 
that increasing g (making the phase function more 
forward throwing) reduces the scattered light levels 
and concentrates the scattered light towards "halos" 
around the UV sources in their simulation. Bianchi, 
Ferrara, & Giovanardi (1996) also pointed out that 
the Heyney-Greenstein phase function may underes- 
timate the forward scattering and this would lead to a 
decrease in scattered light levels. Rather than repeat 
similar calculations for dust scattering of Ha, we note 
the effects of changing the albedo and phase function 
and keep the dust parameters fixed throughout our 
simulation. 

2.4. Radiation Transfer 

In our Monte Carlo radiation transfer we now use 
a scheme that includes "forced first scattering" and 
a "peeling off" procedure as we force the photons to- 
wards the observer with the appropriate weights (see 
Appendix) . We have adopted the forced first scatter- 
ing in order to determine the scattered light contri- 
bution from scattering in optically thin, high latitude 
regions. Our technique then determines the relative 
contributions of direct and multiple scattered Ha ra- 
diation from the WIM and H II regions. 

One of the major questions in the study of the 
WIM is the source of its ionization. Current models 
(Miller & Cox 1993; Dove & Shull 1994; Dove et al. 
1999; Elmegreen 1997) suggest that a clumpy ISM 
may allow photons to leak from the midplanc H II 
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regions to the high latitude gas. While this paper is 
not investigating this problem, we are interested in 
determining how far the Ha probes into the Galaxy. 
To determine the distance that photons have traveled, 
we therefore construct three all-sky images. The first 
contains all photons, representing the total Ha inten- 
sity observed at Earth for this model. The two other 
images contain photons that originated from > 3kpc 
and > 6kpc — using the Monte Carlo technique it is 
straightforward to "label" photons according to point 
of origin, distance traveled, number of scatterings, 
etc. Thus we can determine how far we "see" into 
the galaxy model. We have also separated out direct 
and scattered Ha photons for the WIM and H II re- 
gions, thereby enabling us to quantify the role of dust 
scattering. 

3. Model Results 

Figure 2 shows our all-sky simulation for the ISM 
model and illumination described above. The bright 
point sources have been scaled so that they are satu- 
rated at the maximum intensity from the direct WIM 
emission in the simulation. In Fig. 3, 18° wide latitu- 
dinal cuts are shown for different Galactic longitudes. 
Each panel shows four curves: the heavy line is a 
esc \b\ law and the three other curves show the total 
Ha intensity (upper), scattered Ha from the point 
source H II regions (middle), and scattered Ha from 
the WIM (lower). This simulation reproduces many 
of the observed features of the Galactic Ha emission. 
In particular, the domination of H II regions at low 
latitudes, the csc|6| law is closely followed at inter- 
mediate latitudes, and the intensity is lower than the 
esc |6| law at high latitudes. These gross characteris- 
tics are all present in the data described by Reynolds 
(1990). 

Departures from the esc |&| law at high latitudes are 
seen not only in the Ha, but also in H I observations 
(Dickey & Lockman 1990). For the H I, the fact that 
we are located in a local bubble is cited as the expla- 
nation for the departure. The local bubble is believed 
to be elongated perpendicular to the Galactic plane 
(Cox & Reynolds 1987, Snowden 1986), so there is 
a larger fraction of the H I column density "carved 
out" at high latitudes yielding the apparent depar- 
ture from the esc |6| law. The same effect will also be 
present at Ha, though at a smaller level because of 
the larger scaleheight of the WIM. However, in our 
simulations we approximated the local bubble as a 




Fig. 2. — All-sky Aitoff projection of our Ha intensity 
simulation in Galactic coordinates (l,b). The Galac- 
tic center direction is at the center of the figure and 
Orion is on the right below the midplane. The bright 
point sources have been saturated at the peak level of 
the direct WIM intensity and the WIM and scattered 
light is shown on a square root stretch. 

sphere of radius 200pc yet we still found departures 
from the esc \b\ law at high latitudes. The reason for 
the departures in our simulations is the larger con- 
tribution of dust scattered Ha at low latitudes. Fig- 
ure 4 illustrates this by presenting a esc |6| law, and 
three curves comprising different components of the 
Ha intensity. The upper curve comprises the (direct 
plus scattered) WIM emission and the dust scattered 
Ha from the point sources; the middle curve shows 
the (direct plus scattered) WIM Ha only; the lower 
curve shows only the direct WIM emission. This fig- 
ure clearly shows how scattered Ha is larger at low 
latitudes and yields the departures from the esc |6| 
law at high latitudes. Note that for the direct WIM 
emission, the esc |6| law is closely followed at high 
latitudes, but dust obscuration is responsible for the 
emission being smaller than esc \b\ towards the mid- 
plane regions (\b\ < 10°). In retrospect the departures 
are obvious, since the esc |&| law is derived by integrat- 
ing along paths in a plane parallel layer and does not 
account for any scattered light component. 

An important feature of our simulations is that we 
may determine the scattered light contribution of the 
Ha intensity. We find that the scattered light from 
midplane H II regions is typically around 10% of the 
total intensity at high latitudes (Fig. 3). Towards the 
midplane the scattered light level is larger due to the 
proximity of the H II regions. However, our assump- 
tions of a smooth constant density in each grid cell 
likely yields an overestimate of the scattered light for 
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Fig. 3. — 18° wide latitudinal intensity cuts through 
the simulation. Solid curve shows a esc |6| law, the 
three other curves show the total intensity (upper), 
scattered light from point sources (middle), and scat- 
tered light from the WIM (lower). 

b < 10°. In reality dust will not survive in the intense 
radiation environment of the O stars and more Ha 
photons will escape these regions without scattering 
than predicted by our model. Our analysis is primar- 
ily concerned with estimating the scattered light at 
high latitudes and a more accurate determination of 
the scattered light at low latitudes requires a three 
dimensional model for the ISM density, in particu- 
lar in the environs of the O stars in our simulation. 
Reynolds (1990) estimated the scattered light con- 
tribution from Galactic H II regions to be less than 
15% based on [S II] /Ha and [O III] /Ha ratios mea- 
sured along sightlines towards H II regions and along 
sightlincs in the diffuse ISM. Our model appears to 
confirm this estimate at high latitudes. The results 
of our simulation are in agreement with the Ferrara 
et al. (1996) model of NGC 891. In their externally 
viewed simulation, they estimated an Ha scattered 
light contribution of around 10% at 600pc from the 
midplane. 

We have further separated the Ha intensity accord- 
ing to how far from the Sun the photons originated. 
Fig. 5 shows panels of latitudinal cuts for the total 
Ha intensity (upper) , and the intensity of Ha photons 
(direct or scattered) originating at distances greater 
than 3kpc (middle) and 6kpc (lower) from the Sun. It 
is clear that this simulation is dominated by photons 
originating within 3kpc of the Sun. Note the absorp- 
tion lane in the latitudinal cuts of photons from large 
distances. This arises because for sightlines above the 
dense dust in the H I layer the optical depth is low 
and we see the extended diffuse Ha from the WIM. 
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Fig. 4. — Latitudinal cuts through the simulation 
showing the esc |6| law, and three curves comprising 
different components of the Ha intensity. The up- 
per curve comprises the WIM emission (direct plus 
scattered) and the dust scattered Ha from the point 
sources; the middle curve shows the WIM Ha (direct 
plus scattered) only; the lower curve shows only the 
direct WIM emission. 

At low latitudes the dust limits how far we can see. 
It should be possible to see this effect using the kine- 
matic information provided by the WHAM survey. 

The effects of dust clumping will modify our results 
that are for a smooth density distribution. In addition 
to departures from the smooth scattered light cuts we 
have presented, Witt & Gordon (1996) showed that 
clumping has the overall effect of enabling photons 
to penetrate to larger distances than for the same 
dust mass arranged in a smooth geometry. The Witt 
& Gordon analysis considered radiation transfer in a 
clumpy medium surrounding a point source, which is 
appropriate for the transfer of Ha photons from our 
point sources and will likely raise the scattered light 
levels. However, since the diffuse WIM emission will 
itself be clumpy we intend to extend the Witt & Gor- 
don analysis to include extended sources of emission 
appropriate for analysis and modeling of the WHAM 
dataset. 

4. Polarization Prediction 

With scattered light levels around 10% (as esti- 
mated from line ratio measurements and now also 
from our simulations) , we would expect this dust scat- 
tered Ha to be polarized. Therefore, in addition to 
the scattered light simulations shown above, we have 
also calculated the polarization arising from dust scat- 
tering. The Ha emission from the WIM and the point 
sources is assumed to be unpolarized and the polar- 
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Fig. 5. — Latitudinal cuts through the simulation 
showing the total intensity (upper), Ha originating 
from more than 3kpc (middle) and greater than 6kpc 
(lower). 



ization arises from dust scattering. We have assumed 
that the polarization arises solely from scattering and 
we have not considered the polarizing effect due to the 
transmission of radiation through regions containing 
aligned dust grains. Our neglect of this transmission 
polarization is valid for the low column density sight- 
lines at high latitude since the transmission polariza- 
tion scales with optical depth (see models by Jones 
1989; Jones, Klebe, & Dickey 1992). 




Fig. 6. — Simulated polarization map arising from 
dust scattering of the Ha intensity. 

Figure 6 shows an all-sky map showing the pre- 
dicted polarization for our simulation. The polar- 
ization is largest in the midplane regions where the 
amount of scattered light from the point sources is 
highest. The familiar ccntrosymmetric patterns are 
evident indicative of polarization arising from scat- 
tered light from point sources, or in this simulation 
form concentrations of point sources. In general the 
polarization is perpendicular to the midplane, as ex- 



Fig. 7. — 18° wide latitudinal polarization cuts. 

pected for polarization due to scattering. Transmis- 
sion polarization, which is perpendicular to the large 
scale toroidal Galactic magnetic field and therefore 
parallel to the midplane (Mathewson & Ford 1970, 
Zweibel & Heiles 1997) will reduce the midplane po- 
larization levels. The combination of low scattered 
light levels and polarization cancelation results in the 
signal to noise of our simulations being too small to 
allow a detailed prediction of the polarization at high 
latitudes. To increase the signal to noise, in Fig. 7 we 
have taken 18° wide latitudinal cuts (as in the previ- 
ous section) and averaged the Stokes parameters over 
longitude to show the average polarization as a func- 
tion of latitude. The polarization is at a minimum 
close to the midplane, where the unpolarized direct 
emission dominates. The polarization then increases 
at latitudes just above the midplane due to the large 
scattered light levels. For higher latitudes, where our 
simulations arc more accurate, the polarization de- 
creases and is predicted to be in the range 0.1% to 
1%. If these predicted polarization levels could be 
measured they would provide another means to probe 
the distribution of dust, the dust-to-gas ratio, and the 
scattering properties of the dust in the ISM. 

5. Summary 

This paper has presented a simulation of the diffuse 
Ha intensity incorporating direct emission from the 
WIM and the dust scattering of Ha photons from the 
WIM and midplane H II regions. We have assumed 
that the scattering dust and WIM are distributed ac- 
cording to a smooth axisymmetric prescription. The 
H II regions are assumed to be point sources coin- 
cident with the known locations of O stars within 
2.5kpc of the Sun with Ha luminosities proportional 
to the ionizing flux from each source. For regions be- 
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yond 2.5kpc, we have placed point sources randomly 
into spiral arms and the molecular ring, such that the 
total Ha luminosity from H II regions is equal to that 
from the WIM. 

Our simulation appears to reproduce the gross 
characteristics observed in the Ha background, namely 
the esc |fo law at low latitudes and departures from 
this at high latitudes. The departures from the esc \b\ 
law are attributed, at least in part, to a larger dust 
scattered Ha component from lower latitudes. We 
find that dust scattered Ha from H II regions con- 
tributes less than 20% of the total intensity, although 
some sightlines towards the midplane dominated by 
H II regions show higher scattered light levels. Our 
models therefore are consistent with the line ratio 
analysis of Reynolds (1990), who estimated the scat- 
tered light contribution to be less than 15% of the 
total intensity. We have also made a prediction that 
the scattered Ha should be polarized at around the 
0.1% to 1% level at high latitudes. 

Preliminary results from the WHAM survey show 
that the Ha emission is not smooth and there are 
many filamentary structures present (Haffner et al. 
1998, Reynolds et al. 1998). Detailed modeling of 
these data to determine dust to gas ratios, scattering 
phase functions, etc., will require a 3D treatment of 
the ISM density and WIM emission, as well as inclu- 
sion of large scale Galactic rotation to simulate line 
profiles for the direct and scattered line emission. Our 
codes are set up for this and future modeling efforts 
will be directed towards developing the fully three di- 
mensional models demanded by the observations. 

KW acknowledges support from NASA's Long Term 
Space Astrophysics Research Program (NAG5-6039). 
We thank Jon Bjorkman, Barbara Whitney, and Karl 
Gordon for discussions relating to the forced Monte 
Carlo technique. 

Appendix — Radiation Transfer Algorithm 

In our previous Monte Carlo investigations we were 
simulating the transfer of photons within an axisym- 
metric medium and formed images for external view- 
ing. Upon exiting the simulation, scattered light im- 
ages were formed by projecting the photons into (x, y) 
locations on image planes. Such a technique is effi- 
cient for axisymmetric, external viewing situations. 
However, for our current investigation where we are 
wanting to form images viewed from within the den- 



sity grid, we cannot use this technique since photons 
will rarely scatter in a direction that intersects the 
observer's location. We have therefore modified our 
algorithm as follows so that we force weighted pho- 
tons (either direct or scattered) towards the observer, 
thus enabling us to form all-sky images. 

Ha photons are emitted either from the WIM or 
from the point source H II regions. Direct photons 
are emitted in a direction towards the observer and 
have a weight 

HWcct = e- T /4nd 2 (2) 

where r and d are the optical depth and physical dis- 
tance from the point of emission to the observer. 

To calculate the scattered light contribution we 
first choose a random direction of travel from 4ir stera- 
dians and calculate the optical depth, n, from the 
point of emission to the edge of our grid. We then 
force the photon to scatter at an optical depth less 
than n, sampled from 

T=-log[l-£(l-e- T1 )] (3) 

this reproduces the correct probability distribution 
for optical depths and ensures that all photons scat- 
ter at least once. The photons are subsequently scat- 
tered into directions that are randomly sampled from 
the scattering phase function (Eq. 1). A new optical 
depth is then generated from t = — log £ and the scat- 
tering location corresponding to this optical depth is 
determined. If the location is outside our grid, the 
photon is terminated, otherwise a new scattering an- 
gle and optical depth are generated and the photon 
is tracked until it exits the grid. At each scattering 
location we "peel off" and direct towards the observer 
a fraction of the photon's energy. The weight of this 
"peeled off" photon is 

W scatt = a N (l - e-^) C -^HG{9)/d 2 2 (4) 

where N is the number of scatterings the photon has 
undergone up to that point, n is the first optical 
depth, and t-i and o?2 are the optical depth and phys- 
ical distance from the scattering location to the ob- 
server. The function HG{9) weights the photon by 
the scattering phase function, where 9 is the angle 
through which we force the photon to scatter towards 
the observer. The total intensity is the sum of the 
weights of the direct (Eq. 2) and scattered (Eq. 3) 
photons. We have tested this algorithm for exter- 
nal viewing situations against several "regular" Monte 
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Carlo codes that do not employ any forced scattering 
or weighting algorithms and find it to be accurate and 
efficient. 
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